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Abstract 



We give an elementary new method for obtaining rigorous lower bounds on the connective 
constant for self-avoiding walks on the hypercubic lattice Z . The method is based on loop 
erasure and restoration, and does not require exact enumeration data. Our bounds are best 
for high d, and in fact agree with the first four terms of the 1/d expansion for the connective 
constant. The bounds are the best to date for dimensions d > 3, but do not produce good 
results in two dimensions. For d = 3,4,5,6, respectively, our lower bound is within 2.4%, 
0.43%, 0.12%, 0.044% of the value estimated by series extrapolation. 

Key Words: Self-avoiding walk, connective constant, loop erasure, random walk, 1/d ex- 
pansion. 

1 Introduction 

An n-step self-avoiding walk on the hypercubic lattice Z d is a sequence to = (cu(0), . . . , to(n)) 
of points in Z d , with io(i) and u(i + 1) separated by Euclidean distance one, subject to the 
constraint that oo(i) ^ <^(j) for i ^ j. Unless otherwise stated we take cu(0) = 0. The 
self-avoiding walk provides a model of a polymer molecule with excluded volume. Also, its 
equivalence to the N = limit of the N- vector model has made it an important test case in 
the theory of critical phenomena. 

Let c n denote the number of n-step self-avoiding walks. It has been known for almost 
forty years [1] that the limit // = lim n _ ) . 00 c^J 71 exists and is finite and positive. Moreover the 
subadditivity argument showing existence of this limit also shows that // = inf n >i c^J 71 . This 
limit is known as the connective constant, and is the analogue of a critical temperature for 
the N- vector model. Roughly speaking, // measures the number of sites available for the next 
step of a long self-avoiding walk. The connective constant is admittedly of lesser interest 
than the critical exponents, because it is lattice-dependent while the critical exponents are 
universal. Nevertheless, much work has been done in finding rigorous upper and lower bounds 
for the connective constant, principally on two- and three-dimensional lattices. A review of 
work through 1982 is given by Guttmann [2]. 

Recently, two of us [3, 4] proved that the critical exponents for self-avoiding walks in 
dimensions d > 5 take their mean-field values (7 = 1, z/ = 1/2, 77 = 0) . One key ingredient 
in this proof was an accurate numerical lower bound on the connective constant //. Unfor- 
tunately, we were unable to prove such a numerical bound with existing methods; in fact 
the previous methods give estimates which deteriorate as d — > 00 (see Section 6). Therefore, 
we were led to develop a new method for obtaining lower bounds on //, using loop erasure 
and restoration. This method (with the improvements presented here) provides bounds that 
agree with the first four terms of the 1/d expansion; for d > 4 they are within 0.43% of the 
best numerical estimate of //, and even for d = 3 they are within 2.39% (greatly improving 
the best previous lower bounds). We therefore thought it worthwhile to develop these meth- 
ods in detail; that is the goal of the present paper. The method presented here involves a 



2 



conceptual simplification of the methods used in [4], and also leads to better lower bounds. 
Remarkably, even the most elementary of our new methods leads to a better lower bound 
in d = 3 than has been obtained previously, including the enumeration bound 4.352 of [2]. 
Table 1 summarizes our best bounds and compares them to previously obtained bounds on 



d 


previous bound 


this work 


estimate 


upper bound 


2 


2.620 02 a 


2.305 766 


2.638 1585 (10) d 


2.695 76 6 


3 


4.437 33 c 


4.572140 


4.683 907 (22) e 


4.756 6 


4 


6.718 00 c 


6.742 945 


6.772 0(5)' 


6.832 6 


5 


8.821 28 c 


8.828 529 


8.838 6 (8) fl 


8.881 6 


6 


10.871 199 c 


10.874 038 


10.878 8 (9) 9 


10.903 6 



Table 1: Rigorous lower and upper bounds on the hypercubic-lattice connective constant //, 
together with estimates of actual values, for dimensions 2,3,4,5,6. Errors in the last digit(s) 
are shown in parentheses. 

a) Conway and Guttmann [5], b) Aim [6], c) Hara and Slade [4], d) Guttmann and Enting 
[7, 8], e) Guttmann [9], f) Guttmann [10], g) Guttmann [11]. 

The evaluation of our bounds requires some numerical computation, for which we have 
obtained rigorous error estimates. All numerical values reported in this paper are accurate 
up to rounding of the last digit, except for lower bounds on fi which have been truncated so 
as to provide true lower bounds. 

Our methods can be applied to SAWs on any regular lattice. But for simplicity we restrict 
attention here to the hypercubic lattice Z . 

The plan of this paper is as follows: In Section 2 we describe the method of loop erasure 
and restoration, and systematize the lower bounds on fi that can be obtained from it. These 
bounds involve generating functions of random walks with taboo sets, and in Section 3 we 
show how these taboo generating functions can be computed in terms of the massless-free- 
field lattice propagator. In Section 4 we discuss some aspects of the lower bounds on fi 
obtained. In Section 5 we remark on a different method for proving lower bounds on 
based on comparison with the Ising model. In Section 6 we show that our best bounds agree 
with the 1/d expansion for fi through order d~ 2 . In Appendix A we summarize our methods 
for the rigorous numerical calculation of quantities involving the free-field lattice propagator, 
and give 1/d expansions for various simple-random- walk quantities. 

In a separate paper [12], two of us give a rigorous 1/d expansion for fi through order d~ 3 } 
along with a similar expansion for the critical point of nearest-neighbour Bernoulli bond 
percolation. 
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2 Loop erasure and restoration 



2.1 Definitions 

To describe our loop-erasure-and-restoration method, we need to introduce a number of 
generating functions. For this we need several definitions. 

An n-step walk (n > 0) is an ordered sequence u = (cu(0), . . . , of points in Z d 

such that each point is a nearest neighbour of its predecessor, i.e. \co(i) — u(i — 1)| = 1 for 
1 < i < n. We denote the number of steps in a walk u by \u\. The walk u is said to be a 
memory-T walk, if oo(i) ^ <^(j) for all z, j satisfying < \i — j | < r. We denote by T (x, y) the 
union over all n = 0, 1, 2, ... of the set of memory-r n-step walks from cu(0) = x to u(n) = y. 
Thus, ft (x } y) is the set of all walks from x to y } ft 2 (x } y) is the set of all walks having no 
immediate reversals, 4 (x,j/) is the set of all walks having neither immediate reversals nor 
elementary squares, and so on. The elements of fi are called simple (or ordinary) random 
walks. We denote by 0(x, y) = n^ O n (x, y) the set of all self-avoiding walks (of any number 
of steps) which begin at x and end at y. Finally, we denote by 0(x,«) = U ye7L dil(x } y) the 
set of all self-avoiding walks (of any number of steps) which begin at x and end anywhere. 

A walk to = (cu(0), . . . , cu(n)) is said to be a loop if cu(0) = u(n). We write C T (x) = 
Cl T (x,x) for the set of memory-r loops starting and ending at x. Note that such loops are 
allowed to pass through x many times, and that C T (x) includes the zero-step walk u = (x). 
Note also that the memory-r constraint does not apply modn: for example, io(n — 1) is 
permitted to equal cu(l). 

Given a nonnegative real number /3, we define the generating function (or two-point 
function or Green function) for memory-r walks, 

C T (x,y;P)= (2- 1 ) 

Lu£Q T (x,y) 

and for self-avoiding walks, 

G(x,y;P)= J2 /3 M - ( 2 - 2 ) 

Denoting the number of n-step memory-r walks (starting at the origin and ending anywhere) 
by c n?T , we also define the susceptibilities 

oo 

Xr (/3)= £ C T (0,a;/3) = 5XT/? n (2-3) 

x ei d n =° 

and 

oo 

X (/3)= J2 G(0,x;P) = '£c n p n . (2.4) 

xei d n=0 

For f3 > the sums (2.1)-(2.4) are always well-defined, although they will be +oo for 
sufficiently large {3. In fact, since lim^oo c^J 71 = \i and c n > we have 

(1 - /3/i)" 1 < x(P) < 00 for < /3 < /x" 1 (2.5a) 
X(P) = oo for p > pr 1 . (2.5b) 
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Similarly, the same subadditivity argument which implies existence of the limit defining // 
can also be used to prove that 



fi 7 



Jim = inf (2.6) 



for all < r < 00. It follows from (2.6) that 

(1 - /3/i T ) _1 < Xr(/3) < 00 for < f3 < (2.7a) 
Xr(P) = 00 for/^/i; 1 . (2.7b) 

Clearly /i > ^2 > ^4 > • • • > f- A subadditivity argument can be used to show that 
lim T ^ 00 // T = // (see for example Lemma 1.2.3 of [13]). Since c„ j0 = (2d) n and c„ j2 = 
2d(2d — l) n_1 , we have /i = 2 J and /i 2 = 2o? — 1. The value of /i 4 is shown in [14] to be given 
by the unique positive root of the equation 

x 3 - 2(d - l)x 2 - 2(d - l)x - 1 = 0. (2.8) 

Although methods are described in [14] by which in principle fi T can be computed for r > 6, 
in practice these methods are difficult to carry out in general dimensions. 

To compute our lower bounds on fi we will need the numerical values of C T (0, x; /i" 1 ) 
for a finite collection of sites x. For r = this is given by the well-known Fourier integral 
("free-field lattice propagator") 

/d^k ^/ik-x 
7 — u — — , (2-9) 
(2ttY 1 _ 2df3D{k) ' V ; 



where 



_ 1 d 

D ( k ) = lE cos k i, k = (k 1 ,...,k d ) . (2.10) 

8 = 1 



This expression is valid for < fi < At the critical point fi = 1 = the integral (2.9) 
is finite for d > 2 but infinite for d < 2. An effective means for computing the numerical 
values of (2.9) to high precision for d > 2 is discussed in Appendix A. 

To study d < 2 (of course it is d = 2 which is of interest here), we introduce the potential 
kernel [15, 16] 

A (x-J) = Co(0,0;/3)-C o (0,x;/3), (2.11) 

which remains finite in all dimensions, as f3 f Indeed, it is an immediate consequence 
of the Lebesgue dominated convergence theorem that A (x) = lim A (x; f3) is given by the 



absolutely convergent Fourier integral 



J ' W l-D(k) 1 ' 
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In dimension d = 1, an easy calculation yields A (x) = \x\. Remarkably, in dimension d = 2 
this integral can also be performed analytically for any x [15, Section 15]: for example, 
Ao(ei) = 1, Ao(2ei) = 4 — 8/tt and A (ei + e 2 ) = 4/7T. [Here ei,e 2 are the canonical unit 
vectors in Z 2 .] 

For r = 2 the two-point function can be evaluated using the identity 



1-/3 2 
1 + (2d- l)/3 2 ' 



^(0, z; P) = , , /OJ -,^ 9 <?o(0, z; /?/[! + (2rf - l)/3 2 ]) , (2.13) 



which was derived using convolution methods in [17]. At the critical point fi = p 2 1 = 
this reduces to 

2d — 2 

C 2 (0, x; ^) = ^T^O, x- . (2.14) 

Unfortunately we do not know how to compute the numerical values of C T (0, /i" 1 ) for 
r > 4; because of this we primarily restrict attention in what follows to r = 0, 2. 

We remark that unlike the finite-memory case, it is believed that the self-avoiding-walk 
critical two-point function G(0, x; /i -1 ) is finite for all x in all dimensions, including d = 1, 2. 
This has been proven for d > 5 in [3, 4] — and of course it is trivial in d = 1 — but it remains 
unproven in dimensions 2, 3, and 4. It has been known for some time that G(0,x; f3) = 00 
for f3 > /i- 1 [18]. 



2.2 Identities 

We would like now to establish an inequality relating the two-point functions (2.1) and (2.2). 
For this, we recall the following loop-erasure algorithm, which has been studied in detail by 
Lawler [16]. Given a walk u £ O r (x,j/), we can associate to it a (typically shorter) self- 
avoiding walk p £ 0(x,j/) by erasing loops in an appropriate sequence. We begin by finding 
the last time t\ such that u{t\) = 00 (0) = x } and then erase the sites cu(l), cu(2), . . . ,u(ti) 
from cu, producing a walk which we call p^K In other words we have erased the largest 
possible loop at the site x, namely L = (u;(0), . . . ,u(ti)). [If u does not visit x more than 
once, then we can think of having erased a trivial loop.] The walk p^ does not visit x 
more than once, but it may visit the site p^\l) repeatedly. Let t 2 denote the last time that 
P^Kh) = an d erase the sites / o^ 1 )(2), . . . , /0^(t 2 ) as before. Note that the erased loop 

L x = (^(l),...,^ 1 )^)) cannot pass through 00 (0) = x. This procedure gives rise to a 
walk p( 2 \ which does not visit p^ 2 \0) or p^ 2 \l) more than once. We repeat this procedure 
successively for p^ 2 \ p( 3 \ etc., until arriving at a result which is devoid of loops, or in other 
words which is self-avoiding. For each r, this defines a one-to-one mapping from O r (x,j/) 
into the set lZ T (x } y) whose elements are of the form (p, L 0} Z 1? . . . , L n ), where p £ 0(x, y) 
is an n-step self-avoiding walk (for some n) and each Li £ C T (p(i)). We refer to p as the 
self-avoiding backbone of u. 

In fact it is not difficult to see precisely what the image of this mapping is, or in other 
words to see exactly which elements of lZ T (x } y) can be produced by this procedure, at least 
for r = or 2. One way to do so is to try to reverse the procedure, by beginning with an 
element of lZ T (x } y) and associating to it the walk u which is given by first following the 
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steps of L 0} then taking the first step of p } then taking the steps of L\ } then taking the 
second step of p } and so on. For r = 0, each possible simple random walk from x to y can 
be obtained in precisely one way by this procedure, provided that the loop Li attached at 
p(i) does not intersect any of the previous sites p(0), . . . , p(i — 1), for all i = 1, . . . , \p\. For 
r = 2 the situation is similar but slightly more involved: the attached loops must again 
avoid the previous sites as above, but in addition the next-to-last site of the loop Li must 
avoid the next site p(i + 1) of the backbone (except of course for i = \p\ } when this constraint 
is vacuous). For memories r > 4 the situation is more complicated, due to the presence of 
inter-loop restrictions, and we refrain from entering into details. 

The above discussion can be summarized with identities, which are stated below for r = 
and r = 2. To state the identities we first define the generating functions with taboo set A: 

Cf(x,yJ) = E /3 M ( 2 - 15 ) 

w£Sl r (j:,!/):wni=0 

C^(x,y;P) = E /3 M - (2-16) 

"£fir(i,!/): ujnA=0 
^(M-1)7? Z 

In both cases the sum is over walks which avoid the set of sites A; in the latter case we impose 
the additional restriction that the next-to-last site of to is not z. Clearly, C^ ,z (x, y; fi) < 
C^(x,y; (3), and both quantities are decreasing functions of r and of the set A. We also 
define p[0,j) to be the set of sites {p(0), . . . , p(j — 1)}, for j = 1, . . . , and let p[0, 0) be 
the empty set. We can then write the identities 

C (x,y;P) = J2 ^ Pl UC^ J \p(j),p(jyj) (2.17) 

peQ(x,y) 3=0 

C 2 (x,y;P) = £ ^UC^^ipU^pU)^), (2.18) 

peQ(x,y) 3=0 

where the j = \p\ term in this last equation should be interpreted as C^ '"^ (since the site 
p(j + 1) is nonexistent). For higher memories it is less straightforward to write the analogous 
identities, because of the inter-loop constraints; but by dropping those constraints we have 
immediately the inequalities 

\p\ 

Cr(x,y;P) < f^\{C^\p{j),p{jyji) forr>0 (2.19) 

peQ(x,y) 3=0 

Cr(x,y;P) < P M UCfjM +1 \p(j),p(jy,P) forr>2. (2.20) 

peQ(x,y) 3=0 

2.3 Inequalities (first version) 

It is certainly difficult to analyze the right sides of (2.17)— (2.20) exactly, but we wish to do 
something less ambitious: we will obtain upper bounds by relaxing the avoidance constraints 
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on the attached loops. In particular, we can obtain upper bounds by replacing the restriction 
that the loop attached at p(j) avoid all of p[0,j) by the weaker restriction that it avoid only 
the smaller set 

p \j-k,j) = {p(j-k),...,p(j-l)}, (2.21) 

where A; is a (small) fixed nonnegative integer. [For k = 0, p[j,j) = 0; and for k > j we 
omit the nonexistent points p(i) with i < 0.] Applying this to (2.17) / (2.19) leads to the 
inequality 

\p\ 

C T (x,y;P)< £ pMUcF- k j\p(j),p(j);l3) (2.22) 

peQ(x,y) 3=0 

for any r, k > 0. Now the set of sites p[j — k,j) is (for j > k) simply the range of a (k — 1)- 
step self- avoiding walk starting at a nearest neighbour of p(j)' : so we can get a further upper 
bound by taking the maximum over all such sets. Taking into account translation invariance, 
this leads us to define 

A T (k;/3) = maxC^(0,0;/3) , (2.23) 

where the maximum ranges over all ^-element sets A which are the range of a (k — l)-step 
self-avoiding walk starting at a nearest neighbour of the origin. This maximum is obviously 
a nonincreasing function of r and k. Clearly we have 1 

Cr(x,y;P) < /3 H A T (0;/3)A T (l;/3)---A T (A ; -l;/3)A T (A ; ;/3)l''l +1 - fc 

peQ(x,y) 

= a Ttk (P)G(x,y;PA T (k;P)), (2.24) 

where 



«T,fc(/3) 

Summing over y £ Z d , this gives 



l-k 



i=0 



A T (k; P) . (2.25) 



XAP) < x(PMk-J)). (2-26) 



&T,k(P) 

This is our fundamental loop- erasure- and-restoration inequality , in its simplest form. 
We shall argue below that 

lim = oo (2.27) 

for r = 0, 2, all k > and all d > 0. We expect that (2.27) is true also for all r > 4, but 
we have not proved this, nor shall we make use of it. Given (2.27), it is now easy to obtain 
a lower bound on p, in the following way. By (2.27) and (2.26), x T ( / u~ 1 A T (k; p~ x )) = +oo. 
But we know from (2.5a) that x( x ) < 00 f° r x < ^/l 1 - We thus conclude that 

A T (k\a~ r ) 1 

y ,Ht ' > - , 2.28 
p T p 



1 The terms with \p\ < k — 1 here are being overcounted on the right-hand side, to simplify the form of 
the inequality. 
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or in other words 

H > ^ for r = 0,2, k > 0, d > . (2.29) 

The bound (2.29) would also follow for r > 4, if (2.27) were proven for such r. 

The claim (2.27) follows immediately from (2.7a) if ct^/^/i" 1 ) is finite. And it follows 
easily from (2.25) that ct^/^/i" 1 ) is finite precisely when A T (0; / u~ 1 ) = C T (0,0; /W" 1 ) is finite 
[recall that A T (z;/3) < A T (0; /3) and A T (k; f3) > 1]. For r = or 2 the finiteness of C T (0, 0; /x^ 1 ) 
in dimensions d > 2 follows immediately from (2.9) and (2.14). This proves the claim (2.27) 
for d > 2. For d < 2, consider first r = 0. We have trivially that A (&; /3) > 1. Also, it is 
well-known (see Section 3.1) that A (l; /3) = Co^(0,0; f3) is uniformly bounded for f3 < 
in all dimensions; hence the same is true for A (i; ^) for all i > 1. Finally, for z = we have 

r const x (1 - 2J/3)-( 2 " d )/ 2 for d < 2 
A o (0; /3) = C o (0, 0; /3) ~ j con st x log(l - 2J/3) for J = 2 (2.30) 

^ const for d > 2 

as /3 | Hence in any dimension d > we see from (2.7a) and (2.30) that A o (0; f3) diverges 
more slowly than Xo(fl) as T ji', therefore, lim j_ xo(fl) / &o,k(fl) = °°- The same argument 

PI 2d 

can be used for r = 2, using also (2.13). 

We remark that another possible approach to bounds on fi of this type, which we will 
not pursue further, would be to substitute some f3 < /i" 1 in (2.26), rather than f3 = /i" 1 , and 
then to optimize over f3. This surely gives an improvement when d = 2 and k = 0, but it 
requires some //-dependent a priori upper bound on x(@)i an( l the only available such bound 
is the very weak Hammersley- Welsh bound [19, 20] (for which explicit constants would be 
required). 

To see explicitly what the bounds (2.29) entail, let us consider first the case r = 0, for 
which /i = 2d. Taking k = and k = 1 gives A o (0; ^) = Co(0,0;^) and A (l; ^) = 
Cq 6 ^(0,0; ^) (by symmetry, where e is any neighbour of the origin), so (2.29) becomes 

" s i' = °i < 2 - 31 > 

2d 

u > —r^ [k = 1] (2.32) 

" ^o {e} (0,0;^) 

For a discussion of the relation between (2.31) and some previously known results, see Sec- 
tion 5. The denominator on the right side of (2.31) is infinite in dimension d < 2, but this 
defect is remedied in (2.32): in Section 3.1 we shall prove the identity 

C o {e} (0,0;^) = 2 - 1 < 2 (2.33) 



C o (0,0; 



2d' 



in all dimensions. This already improves the method used in [4], where an infinite denom- 
inator was always encountered for d < 2. In the trivial case d = 1, (2.32) gives the exact 
answer // > 1. For d = 2, (2.32) gives the rather poor bound // > 2. For J = 3, (2.32) already 
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does better than all previously known bounds, yielding p > 4.475817... . In Section 3.1 
we show how to carry out the computations at least in principle for arbitrary values of k, 
and we give explicit numerical results for k = 0, 1, 2, 3, 4 in dimensions d = 2, 3, 4, 5, 6. The 
resulting lower bounds on p are tabulated in Table 2. In Section 6.3 we study the behaviour 
of these bounds as d — > oo. 

Next let us evaluate the bound (2.29) with memory r = 2. Here p 2 = 2 J — 1. For k = 0, 
we conclude from (2.14) that 

2d-l _ (2d-l) 2 1 
11 ~ " 2d -2 C o (0,0;^) • ^ 

This bound is nontrivial for d > 2, and is a factor 

(2d-l) 2 1 

2d(2d-2) 2d(2d-2) [ ] 

better than the corresponding bound (2.31) based on r = 0. For k = 1 the situation is less 
simple: we will show in Section 3.2 that Cj e ^(0,0; ^Zi) can be obtained by solving a linear 
system of three equations in three unknowns. In Section 3.2 we show in fact how to carry out 
the computations at least in principle for arbitrary values of k, and we give explicit numerical 
results for k = 0, 1, 2, 3, 4 in dimensions d = 2, 3, 4, 5, 6. The resulting lower bounds on p are 
tabulated in Table 2. In Section 6.3 we study the behaviour of these bounds as d — > oo. 



2.4 Inequalities (second version) 

The foregoing inequalities are based on constraining the attached loops Li to avoid the 
preceding k sites of the backbone p. For any memory r > 2, we can improve these results by 
using (2.18) / (2.20), i.e. by taking into account the further constraint that the next-to-last 
site of the loop avoid the next site of the backbone. The analysis given previously can be 
repeated almost verbatim in this case. We introduce 

A T (tJ) = maxC^' e (0,0;/3) , (2.36) 

A,e 

where the maximum ranges over all ^-element sets A which are the range of a (k — 1)- 
step self-avoiding walk starting at a nearest neighbour of the origin, and over all nearest 
neighbours e of the origin satisfying e ^ A. The analogue of (2.24) becomes 2 

Cr(x,y;/3) < /3 |p| A T (0;/3)---A T (A ; -l;/3)A T (A ; ;/3)l''l- fc A T (A ; ;/3) 

peQ(x,y) 

= a T , k (f3)G(x } yJA T (k;f3)) } (2.37) 

where 

"tM = (li A T (£; f3)~ k . (2.38) 

2 Here we are overcounting by neglecting the avoidance of the next backbone site, for the first k sites of 
the backbone, to allow for a unified treatment of \p\ < k — 1 and \p\ > k. 
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Arguing as before, we have 

// > ^ — — for r = 2, k > 0, d > (2.39) 

(and again we expect, but have not proved, that this bound holds for all r > 2). 

We will apply (2.39) with r = 2, which requires evaluation of C^' e (0, 0; /i^" 1 ). In Sec- 
tion 3.3 we will prove the identity 

Cf e (0, 0; 0) = Y^iCiiO, 0; /3) - /3C 2 A (0, e; /?)], (2.40) 

from which the bound (2.39) can be computed once we know the values of C^(0, 0; p^ 1 ) and 
Cf(0,e;p?). 

The use of (2.40) involves two improvements on our earlier method [4]. The first im- 
provement is conceptual, in that the derivation here of (2.40) is simpler and avoids the 
combinatorial niceties encountered in [4]. The second is that our earlier method used only 
\A\ = 1, and moreover constrained loops to avoid the previous backbone site only with their 
first and next-to-last sites, rather than entirely as in the methods of this paper. 

Our method for numerically evaluating the denominators A (k; /3), A 2 (k; f3) and A 2 (k; f3) 
which appear in the lower bounds on p is described in Section 3. The numerical bounds 
resulting from (2.29) with r = 0, 2 and (2.39) with r = 2 are given in Table 2 as a function 
of the parameters r and k. We have restricted attention to r = and r = 2 due to our 
inability to compute the numerical values of A T (k; f3) or A T (k; f3) for higher memories. 

In Section 2.5 we will show how to improve on the bounds obtained so far. 



2.5 Optimized bounds 

Our method thus far has been based on taking the maximum over possible geometries for 
the incoming walk (i.e. the set A) in (2.23) or (2.36). This procedure is costly, because 
typically we expect that the incoming two steps of the backbone will be bent rather than 
straight, but the maximum for the small values of \A\ we are using corresponds to a straight 
backbone. In fact, given that there are on average about p possible steps for a self- avoiding 
walk, the proportion of straight to bent steps should be roughly one to p — 1. (This is not 
exactly right, because the straight and bent steps have different probabilities for respecting 
self- avoidance; but it does indicate the expected order of magnitude.) In this section we show 
a way of partially taking this into account to obtain an improved bound. For concreteness, 
we consider the case of memory r = 0, with k = \A\ = 2. 
Summing the identity (2.17) over y £ Z , we have 

Xo(/3) = E ^ II C^\p(j), p(j); P) . (2.41) 

p£fi(0,.) 3=0 

Relaxing the avoidance constraints on the attached loops, to k = 2, we obtain 

Xo(/3) < £ /3 |p| II c^\ P (j), pU); P) ■ (2-42) 

p£fi(0,.) 3=0 
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( rr U\ 

{T,k) 


a = Z 


A Q 

a = o 


A A 

a = 4 


A K 

a = o 


A P. 

a = D 


(0,0) 





3.956 775 


6.454 386 


8.648 213 


10.743 414 


(0,1) 


2 


4.475 817 


6.704 650 


8.809186 


10.862 525 


(0,2) 


2.161367 


4.492 416 


6.706 478 


8.809 464 


10.862 584 


(0,3) 


2.234 696 


4.499126 


6.706 931 


8.809 501 


10.862 588 


(0,4) 


2.275 515 


4.502 565 


6.707 091 


8.809 509 


10.862 588 


(2,0) 





4.121 641 


6.588 853 


8.756 316 


10.832 942 


(2,1) 


1.712626 


4.441266 


6.696 516 


8.806 308 


10.861210 


(2,2) 


1.917318 


4.457137 


6.698 205 


8.806 565 


10.861265 


(2,3) 


2.019 782 


4.463 834 


6.698 632 


8.806 599 


10.861268 


(2,4) 


2.079 398 


4.467 350 


6.698 783 


8.806 607 


10.861269 




n 




U.UO 1 <JO<J 


o. i ou uoy 


iu.o4:U uuy 


(2,1) 


1.976 372 


4.539 419 


6.737460 


8.827159 


10.873 577 


(2,2) 


2.153 350 


4.552 467 


6.738 907 


8.827 387 


10.873 627 


(2,3) 


2.239 265 


4.557 994 


6.739 273 


8.827417 


10.873 630 


(2,4) 


2.286 245 


4.560 903 


6.739 404 


8.827424 


10.873 631 


(0, 2) 0?rf 


2.195 201 


4.518 652 


6.715 924 


8.813 103 


10.864 240 


(0, 3) op t 


2.267128 


4.526 286 


6.716 713 


8.813 204 


10.864 257 


(OA) opt 


2.305 766 


4.530 282 


6.716 982 


8.813 224 


10.864 259 


(2, 'ijopt 


1.936810 


4.476 092 


6.704487 


8.808 707 


10.862112 


(2, 3) pt 


2.038 216 


4.483 773 


6.705 232 


8.808 801 


10.862128 


(2, 4:)opt 


2.092 741 


4.487 869 


6.705 488 


8.808 820 


10.862130 


(2, 'ijopt 


2.165 878 


4.562 269 


6.742 085 


8.828 430 


10.874 022 


(2, 3) pt 


2.248 707 


4.568 677 


6.742 721 


8.828 512 


10.874 036 


(2, 4) pj 


2.290 302 


4.572140 


6.742 945 


8.828 529 


10.874 038 



Table 2: The lower bounds on /i of (2.29) for r = 0, 2 and of (2.39) for r = 2, for k = 1, 2, 3, 4, 
and the optimized bounds using the method of Section 2.5, for dimensions d = 2,3,4,5,6. 
The numerical values have been truncated to give rigorous lower bounds. The lower line of 
the table provides the best lower bounds, except for d = 2. Some discussion of the table 
entries is given in Section 4. 
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Our previous approach was to replace the quantities by the corresponding upper 

bounds based on using the worst 2-element set A in place of p[j — 2, j) [cf. (2.23)/(2.24)]. 
Now let us try instead to distinguish the two possible cases, 

A = straight (i.e. congruent to {ei,2ei}) 
A = bent (i.e. congruent to {ei,ei + e 2 }). 

Correspondingly, let us define 

a = C o {ei ' 2ei} (0,0;^) (2.43a) 
b ee C o {ei ' ei+e2} (0,0;^). (2.43b) 



2d 

2d 

We only consider the case of 



< b < a < 2d } (2.44) 
which is shown to be valid for our applications by numerical computation. It follows that 3 

Xo(/3) < £ /3 W Ao(0;/3)Ao(l;/3)a^)6l"l- 1 -^) , (2.45) 
P efi(o,.) 

where £ denotes the number of straight vertices among p(l), . . . , p(n — 1) [here n = \p\]. 
Defining 

c n (£) = £ 0(0, •): \p\ = n, ^(straight vertices in p) = £} , (2.46) 

we can rewrite this as 

oo n — 1 

Xo(/3) < Ao(0;/3)A o (l;/3)EE c 4^)/3^" 1 "' 

n=0 1=0 

oo 

ee A o (0;(3)A o (l;(3)J2l3 n S(n;a,b) (2.47) 

n=0 

where 

n-l 

S(n;a,b) ee ]T c n (t)a l b n - x - 1 . (2.48) 
i=o 

As in the preceding subsections, we will let f3 f and use the fact that 

This together with (2.47) implies that 

limsup5'(n;a,6) 1/n > 2d. (2.50) 



3 The terms with \p\ = 0,1 are here being overcounted on the right-hand side. 
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Therefore, if we can get upper bounds on S(n; a, b) in terms of the {c n } } then we will be able 
to deduce lower bounds on //. 

Example 1. Since a > b > 0, we have trivially 

S(n;a,b) < c n a n ~ x (2.51) 

and hence 

lim sup S(n;a,b) 1/n < jia . (2.52) 

Combining this with (2.50), we obtain // > 2d /a. This is, of course, our old bound (2.29) 
with r = 0, k = 2, based on using the worst set A. 

Example 2. We claim that 

c n (£) < 2di^~ l ^jl t {2d-2) n - 1 - t . (2.53) 

Indeed, there are ways of distributing I straight vertices among the n — 1 internal 

vertices of an n-step walk; and at each straight (resp. bent) vertex there are at most 1 (resp. 
2o? — 2) choices for the next step. The first step of the walk has, of course, 2o? choices. [We 
remark that summing (2.53) over i gives c n < 2d(2d — l) n_1 , which is the trivial bound on 
c n in terms of memory-2 walks.] Inserting this bound into (2.48) and performing the sum, 
we conclude that 

S(n;a,b) < 2d[a + (2d - 2)b] n ~ 1 (2.54) 

and hence 

limsup5'(n;a,6) 1/n < a + (2d - 2)b . (2.55) 
Combining this with (2.50), we learn that 

a + (2d-2)b>2d (2.56) 

— an interesting fact, but one which unfortunately teaches us nothing about //. 

Our approach now will be to combine , in an optimal way, these two ways of bounding 
S(n; a, b). We shall use the first bound for £ < Xn and the second bound for £ > An, with a 
suitable choice of A. That is, we shall split S = S\ + S2, with 

[\n\ 

S^n-a.b) = c ^V hn ~ X ~ l (2.57a) 

n-l 

S 2 (n;a,b) = ]T c^aV" 1 ^. (2.57b) 

Since a > b > 0, we have trivially 

Si(n;a,6) < c^^b 11 - 1 -^ (2.58) 



14 



and hence 

limsup5'i(n;a,6) 1/n < //a A & 1_A . (2.59) 

n—^co 

On the other hand, from (2.53) we have 

S 2 (n;a,b) < 2d ]T ( U ~ l \ a 1 [(2d - 2)b] n - 1 ~ t . (2.60) 

£=[\n\+l V / 

Using Stirling's formula, we find that the summand achieves its maximum at 

£ . a , 
- k, A = -7 2.61 

n a + (2d-2)b y ' 

and is decreasing thereafter. So let us take A £ [A , 1]; then easy estimates show that 

limsup£ 2 (n;a,&)^ < ^ ' ^ 

Combining (2.50) with (2.59) and (2.62), we conclude that 



max 



jia b 



A 



(2d-2)b 



l-A' 



XJ \ l-A 



> 2d (2.63) 



for all A £ [A , 1]. We will choose A £ [A , 1] such that the second expression is just barely less 
than 2d; then we can conclude that the first expression is > 2d. Now, the second expression 
equals a + (2d — 2)b > 2J [the inequality due to (2.56)] when A = A , and equals a < 2d 
[the inequality due to our assumption (2.44)] when A = 1, and is a continuous and strictly 
decreasing function of A in the interval [A , 1]. Therefore, there is a unique A* £ [A , 1] such 
that 



a\ A * { (2d - 2)b\ 1 ~ X * 



.A*/ \ l-A* 
We then have the bound 



2d . (2.64) 



2d (2d-2\ 1 A * / 1 x A * 



" * = (—) IfJ • (2 - 65 » 

The above inequality can be used to improve on our previous bounds. For \A\ = 2, it can 
be applied directly. To apply it when \A\ > 2, we classify A according to its first two steps, 
and distinguish the two possibilities by defining now [we write A = (ai, a 2 , 03, ...)] 



max C o A (0,0;^ 

A:ai=ei ,a 2 = 2ei U V 2(1 



max C£(0,0;±) (2.66a) 

A: a 1 =e 1 ,a 2 = e 1 +e 2 



Then (2.45) still holds for these a and 6, and we can proceed in exactly the same way as for 
the \A\ = 2 case. The bounds for memory-2 and memory-2 can be optimized in the same 
way. The resulting bounds are given in Table 2. 
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3 Evaluation of loop generating functions 



3.1 Simple random walk 

In this section we shall reduce the computation of A (k; f3) to the evaluation of the simple- 
random-walk two-point function Co(0,x; f3) at a suitable finite set of sites x. In particular, 
for computing the numerical bounds in Table 2 we need only the critical case f3 = 1/2 d. 
Methods for obtaining the values of the critical two-point function are discussed in detail in 
Appendix B of [4], and a brief summary is given in Appendix A below. Of course the critical 
simple-random- walk two-point function is infinite in dimensions d = 1,2; we briefly indicate 
the modifications needed to treat d < 2 at the end of this section. 

We begin by deriving a recursion relation which relates the generating functions C^(j/, x; f3) 
and CQ U ^(y, x; f3), where b is a single site and A is a finite set of sites which does not contain 
b. Applying inclusion-exclusion gives 

Cf u{ %,x;P) = C£{y,x;j)) - £ /?W . (3.1) 

uj 3 b,uj tl A = 

Given a walk u contributing to the sum on the right side, we break the walk into two pieces 
at its first visit to b. The generating function for walks which go from y to 6, which hit b only 
once (namely at their last step), and which avoid the set A, is equal to C^(j/, b; f3)/C A (b, b; f3). 
Therefore 

r ,Au{b} { r A ( a\ C^ (y } b; f3) C^jb, x; f3) 

C xs {y } x;f3) = C (y } x;f3) fe . (3.2) 

When A is the empty set, the right side can be computed in terms of the ordinary two-point 
function Co, and then by iteration we can compute C^(j/, x; fi) for any finite set A. 

An amusing special case is y = x = and A = {e}, where e is a nearest neighbour of the 
origin. Using the identity 

C o (0,0;/3) = 1 + 2J/3C o (0,e;/3) (3.3) 

together with (3.2), we find 

This is finite as f3 f ^ in all dimensions d > [since Co(0, 0; f3) diverges as f3 f ^ more 
slowly than (1 — 2d/3)~ 1 ] and yields 

In particular, we conclude that A (i; /3) < A (l; /3) = Cq 6 ^(0,0; f3) < 2 for all i > 1 and all 
ft < in all dimensions. 
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The computation of A (k; f3) for any chosen f3 is now reduced to a finite amount of labour. 
We simply list all the allowable sets A of the given cardinality k, exploiting the obvious lattice 
symmetries, and then compute C^(0,0;/3) for each such A by iterating (3.2). For \A\ = 1, 
we have only one choice, A = {e}. For \A\ = 2, we compute the maximum over the two 
choices A = {ei,2ei} and A = {ei,ei + e 2 }, where e 4 - are the canonical unit vectors. And 
so forth. Values for C^(0,0; f3) are tabulated in Table 5 in the Appendix, and the resulting 
lower bounds on // are given in Table 2. 

In dimension d > 2, we can perform these computations directly at the critical point 
ft = 2d : a ^ quantities are finite. However, in dimension d < 2 the two-point function 
Co(0, x; ^) is infinite for all x. Nevertheless we can show that (0, x; ^) is finite whenever 
A ^ 0, by working first at f3 < and then letting f3 | Note first that by the monotone 
convergence theorem we have 

C^x;j- d ) = limC^x-J). (3.6) 

Now the trick is to rewrite the case A = of (3.2), which is certainly valid for f3 < in 
terms of a quantity that (unlike Co) stays finite as f3 f Such a quantity is the potential 
kernel A (x; f3) = Co(0,0;/3) — Co(0,x;/3) discussed in Section 2.1. Rewriting (3.2) with 
A = in terms of the potential kernel gives 

Ci%,x;fl) = -A (x-y; P) + A (b-y; P)+A (x-b; /?)- A ° (& ~ V '] f } ^ ~ b ' J) } (3.7) 

which in the limit f3 | ^ in dimension d < 2 reduces to the result of Proposition 11.6 of [15]: 

C {6} (j/, x; ±) = -A (x -y) + A (b -y) + A (x - b) . (3.8) 

In particular this shows that C^(j/, x; ^) < oo whenever A ^ 0. Thus, once we have handled 

the step from Co to Cq^ using (3.8), we can then advance to larger values of \A\ using (3.2) 
directly at f3 = just as in higher dimensions. (An alternate method of computation 
is given in Theorem 14.11 of [15], but we have found the recursion easier to implement 
numerically.) 

3.2 Memory- 2 walk 

The computation of our memory-2 lower bounds on // has now been reduced to the evaluation 
of the generating functions Cf^O;//^ 1 ) and Cf^e;//^ -1 ) for finitely many finite sets A, 
just as was described in Section 3.1. The basic idea for the evaluation of these quantities 
is the same here as for simple random walk, but the recursion relation requires more care. 
Suppose that b^A. If b = y } we just have 

C? u{6} (y,a;/3) = 0. (3.9) 
lib ^ y, we begin as for simple random walk by writing 

C? u{b} (y,x-J) = C 2 (y,x-J)- ft H - (3-10) 

^ 6 ^2 (y, x) 
uj 3 b,uj tl A = 
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To deal with the sum on the right side, we again want to cut the walk at the first time it hits 
6, but now the two pieces of the walk are no longer independent because of the memory-2 
constraint. The sum on the right side is equal to 

G Q2(y,b),uj 2 G 2 (6,a:) 
uj\ , LU2 n A = 

= £ E /3 KI I[^&] £ ^IM1)^6 + /], 

/:|/|=1 ^en 2 (2/,6+/) c 2 efi 2 (&,:r) 

wJnA = w 2 nA = 

(3.11) 

where u\ o cu 2 denotes the concatenation of the two walks. The sum over uj' x exactly gives 

ct u{b} (y,b + f-j). 

For the sum over cu 2 , we use repeated inclusion-exclusion to remove the constraint that 
cu 2 (l) 7^ b + /, to obtain 

J2 ^i[u 2 (i) ^ b + f] = ci(b,x-j)-pci(b + f, x -j) 

+/3 2 C 2 4 (6, x; /3) - p 3 C?(b + f,x;P) + ... 



w 2 £ £2 2 (&, a;) 
u; 2 n A = 



:[Cf(b,x;P)-pCf(b + f,x;P)]. 



1 - (3 2 

(3.12) 



The result is 

Ci u{h \y,x;j)) = Cf(y,x;P) 



-V^T2 E C' 2 Au{ %,6 + /;/3)[C' 2 4 (6,x;/3)-/3C' 2 4 (6 + /,x;/3)]. 

P /:|/l = l 

(3.13) 

The above unfortunately is not a closed-form expression for the left side, because of the 
occurrence of similar quantities on the right side. However, for a fixed j/, (3.13) provides a 
systems of 2d linear equations for 2d unknowns, namely {C^^ (y , b + /; ft)}\f\=i- Once we 
know these 2d quantities by solving the equations, then everything on the right hand side of 
(3.13) is known, and C^ u (y, x; f3) can then be computed for general x. Fortunately, the 
number of unknowns is often reduced by symmetry. For example, for y = 0, A = and 
b = ei, there are three inequivalent values of b + /, namely 0, 2ei and e\ + e 2 . Thus we can 
first solve the system of three equations in three unknowns which results by taking x to be 
each of the three values 0, 2e\ and e\ + e 2 , and once this system has been solved we can then 
compute the left side directly for any other value of x. Concretely, for d = 3, we get 

C 2 {ei} (0,0;|) « 1.125805 (3.14a) 

C 2 {ei} (0,2 ei ;|) « 0.074441 (3.14b) 

C 2 {ei} (0, ei + e 2 ;i) « 0.138449 (3.14c) 
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For \A\ > 1 we can proceed similarly; results are given in Table 6 in the Appendix. 

The right side of (3.13) can be used directly at the critical point in more than two 
dimensions, but it contains infinities for d = 2 when A is the empty set. In two dimensions 
a limiting argument is first used to deal with A = 0, and then for \A\ > 1 the recursion can 
be used directly as above. Surprisingly, the resulting bounds are worse than those obtained 
using simple random walk, and so we do not give the details of the limiting argument, but 
instead give only the result: for b ^ j/, 

C?W;D = I + |[A (6-y)-A (x-y)] + A (x - b) 

E Ci b} (y } b + f;l)A (x-b-f). (3.15) 

/:|/l = l 

Results are discussed in Section 4.2. 



3.3 Memory-2 loop with taboo on penultimate site 



In this section, we derive the identity (2.40), which shows how the evaluation of C^' e (0, 0; /3) 
can be reduced to the evaluation of C<f (0, 0; f3) and C<f (0, e; f3). By definition, 

Cf e (0, 0; 0) = Ci(0, 0; /3) - £ ^ . (3.16) 

w e n 2 (o,o) 

w n A = 0, - 1) = e 

By applying the inclusion-exclusion relation repeatedly as in (3.12), we obtain 

£ P H = 73^^(0, e; V ~ ^(0, 0; /?)]• (3.17) 

wCi A = 0, - 1) = e 

Inserting this result into (3.16) gives 

Cf e (0, 0; 0) = Y^iCiiO, 0; 0) - /3C 2 A (0, e; /3)], (3.18) 

and hence the bound (2.39) can be computed once we know the values of C^(0, 0; /i^ 1 ) and 
C^(0, e; /U^ 1 ). For \A\ > 1, (3.18) can be used directly at the critical point in all dimensions. 



4 Discussion of results 
4.1 Three or more dimensions 

From Table 2 it can be seen that for d > 3 the r = 2 bound does better than r = when 
k = 0, as was already shown at the end of Section 2.3. Perhaps surprisingly, for k > 1 this 
situation is reversed: higher memory means a more "sophisticated" bound, but it does not 
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necessarily mean a better bound! In any case, the bounds with memory-2 do better than the 
corresponding memory-0 bounds. 

In all cases with memory-0 or memory-2 in Table 2, the set A giving the maximum ol 
Cq(0, 0; ^) is the straight line segment, as can be expected intuitively for small \A\. However 
for large \A\ in general it is not to be expected that the straight segment is optimal, and an 
example of an optimal A which is not straight is given below for d = 2. (We thank Greg 
Lawler and Alain Sznitman for discussions on the non-optimality of straight A for large \A\.) 
In Table 2 it is also the case that straight A gives the maximum for the memory-2 bounds, 
although a priori there seems to be no compelling reason to expect this to be the case: in 
(3.18) straight A will give the maximum individually for each of C^(0, 0; f3) and C^(0, e; 
but possibly not for the difference. 

4.2 Two dimensions 

Unfortunately, for d = 2 our methods do not produce good bounds on //, as can be seen 
from Table 2. For example, with \A\ = 4 and r = we have only // > 2.275515. To 
get a measure of the inherent limitation of the method in two dimensions, and in view 
of the relatively slow convergence of the bounds as k increases, we computed the value of 
2d/CQ k (0, 0; ^) with Lj~ the collinear set of k sites joining ( — 1,0) to ( — k,0). For large k 
we do not expect that C^(0,0; ^) is maximized over sets A with \A\ = k by Co*(0,0; ^), 
and in fact already for k = 5 and A = {( — 1, 0), (— 2, 0), (— 3, 0), (— 4, 0), (— 4, 1)} we have 
C o A (0,0;^) = 1.739044 > 1.738131 = C o i5 (0, 0; ^). Nevertheless 2d/Co" (0, 0; ^) does 
provide an upper bound on the best possible r = (0, k) lower bound on //. For A; = 60 we 
found 2d/Co 60 (0, 0; ^) = 2.404210, and by extrapolating the sequence for k < 60 to k — > oo, 
we found a limiting value less than 2.42. Substantially better lower bounds on // can be 
obtained by other methods [2, 5]. 

The situation for r = 2 is only slightly better. Solving the system (3.15) of three equations 
in three unknowns gives 

C 2 {ei} (0,0;|) « 1.751695 (4.19a) 
C 2 {ei} (0,2 ei ;|) « 0.649131 (4.19b) 
C 2 {ei} (0, ei + e 2 ;|) « 0.799587 (4.19c) 

This gives the very weak bound // > 3/1.751695 > 1.712626. By imposing the additional 
restriction that there should be no direct returns at the end of the loop [the memory-2 bound 
of inequality (2.39)], this is improved to // > 1.976372, which is comparable to the result 
obtained with memory-0. 

Again for memory-2 and memory-2 we computed the bounds arising from A = for 
large k, to obtain an idea of the inherent limitation of the method. For memory-2 and A; = 58 
we found a bound of 2.268661, while for memory-2 and A; = 58 we found a bound of 2.443124. 
Only marginally higher values result when these are extrapolated to k — > oo. Since straight 
A in general will not be optimal, the best possible bound on // using memory-2 or memory-2 
may in fact do worse than these values. In particular, for the memory-2 bounds straight 
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A already fails to provide a maximum when k = 4, where the optimal A is the bent set 
A= {(-1,0), (-2,0), (-3,0), (-3,1)}. 

Because of the poor results for d = 2 we have not done a rigorous error analysis of these 
d = 2 computations, but we do expect that they are correct to the given accuracy. 



5 Comparison to the Ising model and the infrared 
bound 

Our simplest (and worst) bound 

2d . . 

" " cw^i) (5 - 1) 

is in fact not new, as it is an immediate consequence of earlier results of Fisher [21] and 
Frohlich, Simon and Spencer [22]. Also, this bound was proven by Lawler [23] for dimensions 
d > 4, via a different perspective on loop erasure. 

Fisher's result is that the two-point spin correlation function (cr x cr y )j of the nearest- 
neighbour ferromagnetic Ising model with inverse temperature J (and zero external field) 
obeys the inequality 

((T x (7y) J < G(x,y; t&nh. J), (5.2) 

from which he concluded that 

jl > COth J c j s i n g , (5.3) 

where J c ,ising is the critical inverse temperature of the Ising model. The infrared bound of 
[22] gives 



Combining these two inequalities yields 



'Co(0,0;4 



H > coth ^ " v ^ 2dJ J . (5.5) 

This bound is an improvement on (5.1), and yields (5.1) when combined with the inequality 
cothx > x~ x (for positive x). 

Table 3 gives the value of coth J c ,ising , using numerical estimates of J c ,ising , and compares 
it with the bounds (5.1) and (5.5) as well as with the best rigorous lower bounds on //. Values 
of Jcjsing are taken from the exact solution J c ,ising = \ log(l + y/2) for d = 2, from the Monte 
Carlo study [24] for d = 3, and from the series extrapolation results of [25] for d = 4 and [11] 
for d = 5,6. If we accept the (nonrigorous for d > 3) numerical values of J c ,ising } then our 
lower bounds are better than the best possible bound that can arise from (5.3) for d > 4, 
but not for d = 2, 3. However, for d = 2 the enumeration bound of Conway and Guttmann 
[5] does better than the Ising-model bound. The fact that our best bounds do better in high 
dimensions than (5.3) can be explained by the fact that our best bounds capture more terms 
of the 1/d expansion for // than does the right side of (5.3); see Section 6. 
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d = 2 


d = 3 


d = 4 


d = 5 


d = 6 


J casing 


0.440 687 


0.221659 5(26) 


0.149 663(34) 


0.113917(7) 


0.092 294(7) 


2<*/<7 o (0,0;£) 





3.956 776 


6.454 386 


8.648 214 


10.743 415 


coth[iC (0,0;^)] 


1 


4.040 663 


6.505 948 


8.686 723 


10.774424 


COtll J cjsing 


2.414213 


4.585 07(5) 


6.7315(15) 


8.8162(6) 


10.865 7(9) 


our best bound on // 


2.305 766 


4.572140 


6.742 945 


8.828 529 


10.874 038 


best bound on // 


2.620 02 


4.572140 


6.742 945 


8.828 529 


10.874 038 



Table 3: Comparison of the Ising-model lower bound coth J c ,ising with the best lower bounds 
on fi. The values of J c ,ising are nonrigorous numerical estimates for d = 3,4,5,6; errors in 
the last digit(s) are shown in parentheses. Other quantities are rounded to 6 digits, except 
for the lower bounds in the last two lines, which have been truncated. 



6 High-<i behaviour of the lower bounds 

The lace expansion has been used [12] to give a rigorous proof that as d — > oo, 

The next term has been given in [26, 27], 

1 2 11 62 . 

" = 2J - 1 - 2 7^-(27^-(27^-(27^ + --- (6 - 2a) 

= "-'-h-w-w-w*-"' (6 - 2b) 

but with no rigorous bound on the error. In this section we study the 1/d expansion for our 
(and other people's) lower bounds on //, and show that our best bound agrees with (6.1) up 
to and including the term of order d~ 2 . 

In Section 6.1 we analyze briefly the J-dependence of some older lower bounds on //. 
It turns out that most of these methods have very poor behaviour as d — > oo; this was, 
in fact, the original motivation for us to develop the loop-erasure- and-restoration method. 
In Section 6.2 we comment briefly on some results of Kesten [20], which are based on a 
precursor of our method. In Section 6.3 we analyze the high- d behaviour of the loop-erasure- 
and-restoration bounds, and show that they capture the first few terms of the 1/d expansion 
of fi. An interesting structure emerges from the comparison of how many terms are captured 
for different pairs (t, In Section 6.4 we carry out an analogous analysis for the bounds 
based on comparison to the Ising model. 
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6.1 High-<i behaviour of the older lower-bound methods 

We assume in this subsection that the reader is acquainted with the definition and funda- 
mental properties of bridges and irreducible bridges; see e.g. [2, 13]. 

Let b n (resp. i n ) be the number of n-step bridges (resp. n-step irreducible bridges) starting 
at the origin and ending anywhere. By convention, b = 1 and i = 0. These counts satisfy 
the renewal equation 

n 

b n = S n ,o + ^2ikb n -k ■ (6.3) 
k=l 

We define for f3 > the generating functions 

oo 

B(P) = J2 b nP n (6-4) 

oo 

I(P) = J2^ n - (6.5) 



n=0 



Now < i n < b n < fi n } and hence it follows from (6.3) that 



B ^) = Y^W) for/3 - /i_1 - (6 - 6) 



Also, it is known that 



< (1 - /3/")" 1 < oo for < /3 < n~ x 
B(J3) { -> +oo for £ T /w" 1 (6-7) 

= +oo for f3 > /i -1 

<PH<1 for < ft < fi' 1 
I(P) {=1 for = /i- 1 (6.8) 

= +oo for P > 

Now suppose that we have enumerated Zi, . . . , zjv- (In practice one does this by enumer- 
ating &]v an d solving (6.3) [2].) Then defining 

N 

I N (P) = ^i n P n , (6-9) 

71 = 1 

we have In(P) < I(P) < 1 for /3 < It follows that 

fj, > n* iN = — , (6.10) 

P*,N 

where /J^jv is the unique positive solution of the polynomial equation In(P) = 1- 

How does the bound (6.10) behave as d — > oo for fixed iV? We are unable to give the 
precise asymptotic behaviour, but we can give an upper bound as follows: Trivially we have 
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in < b n < (2d — because the first step of a bridge must be in the +ei direction, and 

for the remaining steps there are at most 2 d — 1 choices each. Therefore, 

N 

In(/3) < J2( 2d - 1 T~ 1 i 3n ■ ( 6 - n ) 

71 = 1 

In particular, the solution /J^jv of In(P*,n) = 1 must be greater than or equal to the solution 

/3**,jv of ^2^ =1 (2d — N = 1. Now, it is easy to see that for each fixed iV, as d — > oo 
we have 

^ N = (2J-l)-( JV - 1 )/ JV [l-0(J- 1 / JV )] , (6.12) 

so that 

= < = (2d - 1)^ N [l + 0(d-^ N )} . (6.13) 

P*,N P**,N 

So this method, for any fixed iV, can never get even the correct leading asymptotic order 
(namely fi ~ d) as d — > oo. 

In practice, this method gives the best currently available bounds for d = 2 [2, 5], but it 
is inferior to the loop-erasure-and-restoration method for d > 3. 

Other early lower bounds on // are those of Rennie [28], 

/i>2J-logJ-3 + y2 + lo g 2-f; . r l0g ! 3 ^ 1 } U] , (6.14) 

j^sJU - lo gU + 1)] 

and Hammersley [29], 

fi > 2d-log(2d- 1) - 1 . (6.15) 

These bounds get the correct first term in the large- d expansion, but numerically they are 
rather poor. 

Some other methods for proving lower bounds on fi in d = 2, 3 are given in [14]; but they 
do not appear to behave well as d — > oo. 



6.2 Kesten's bounds 

A precursor of our method was used by Kesten [20] to prove that 

fi> fi 2r - 0(d~ r ) (6.16) 

for all r > 0. Kesten's method involves loop erasure and restoration at the level of counts 
rather than generating functions. Since trivially fi < fi 2r} (6.16) implies that the 1/d expan- 
sion of fi 2r agrees with the first r + 1 terms of the 1/d expansion for fi. Unfortunately, from 
the point of view of numerical estimates on fi } (6.16) is not very helpful, since it is difficult 
to get good constants in the error term. Nor is it easy to compute the 1/d expansion of fi 2r 
for r > 3, so as to obtain the 1/d expansion for fi beyond three terms. However, using r = 2 
in (6.16) together with (2.8), one obtains 

fi>2d-l-^- + 0(d- 2 ) , (6.17) 
2d 

which gives a bound agreeing with the first three terms of the 1/d expansion for fi } albeit 
without good control of the error term. 
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6.3 1/d expansion for the loop-erasure-and-restoration lower bounds 

We now turn to the computation of the 1/d expansions for some of our lower bounds on //. 
We use the shorthand s = The standard of comparison for all our bounds is the series 
(6.2), 

// = s" 1 - 1 - s - 3s 2 - 16s 3 - 102s 4 + . . . , (6.18) 

which is provably correct through order s 3 and presumably correct also at order s 4 . As 
always, we classify our bounds according to the memory they use (r = 0, 2 or 2) and how far 
back on the backbone they enforce avoidance of the attached loops (k = 0, 1, 2, 3, . . .). We 
introduce the notation ^ T,k ^ to denote the bound obtained by considering (r, &)-quantities. 
Our most basic bound (2.31) is based on r = 0, k = 0; using (A. 17), its 1/d expansion is 

9d 

= m = s- 1 - 1 - 2s - 7s 2 - 35s 3 - 215s 4 + 0(s 5 ) . (6.19) 
^o(0,0; Y d ) 

This gets the first two terms correct, but misses the term of order s. The next simplest 
bound (2.32) is based on r = 0, k = 1; using (A. 17) and (2.33), we obtain 



(o,i) 



2d 



- 1 - 1 - s - 4s 2 - 22s 3 - 143s 4 + 0(s 5 ) . (6.20) 



This gets the first three terms correct, and just barely misses the term of order s 2 . As can 
be seen from Table 8 in the Appendix, taking k = 2 or 3 does not improve on (6.20). 

Next let us consider the bounds based on the memory-2 loop. The simplest such bound 
(2.34) is based on r = 2, k = 0; using (2.14) and (A. 17), it becomes 



(2,0) 



2d-l 



^(0,0;^) 



s" 1 - 1 - s - 6s 2 - 35s 3 - 222s 4 + 0(s 5 ). (6.21) 



This gets the first three terms correct. The bound is further improved if we go to k = 1: 
making use of Table 9, we have 



(2,1) 



2d-l 



s" 1 - 1 - s - 4s 2 - 23s 3 - 150s 4 + 0(s 5 ). (6.22) 



c 2 {e} (o,0;^r) 

The term of order s 2 is improved compared to (6.21), but the coefficient is still not correct. 
This bound based on r = 2, k = 1 is inferior (at order s 3 ) to the simpler bound based on 
r = 0, k = 1, as was the case for the numerical values of Table 2. Using Table 9, it can be 
seen that the expansion for fji 2 ' 2 ^ is identical to (6.22). 

Better memory-2 bounds can be obtained by insisting that the attached loops avoid also 
the next site on the backbone (Section 2.4). The simplest such bound is 



(2,0) 



2d-l 



- 1 - 1 - s - 5s 2 - 29s 3 - 188s 4 + 0(s 5 ), (6.23) 
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where we have used Table 10. This gets the first three terms correct, and does better at order 
s 2 than (6.21), but still has not yet got the correct coefficient —3. This coefficient is however 
captured correctly if we combine the constraints involving the previous and next steps on the 
backbone, i.e. if we go to r = 2, k = 1 . To see this, consider the loop generating functions 
£<{ei}, e i^Q^Q. ^-J— j-) and Cj 61 ^' 62 (0, 0; ^~[), which by symmetry are the only two geometries 
to be considered. Intuitively, the second should be greater than the first, since a walk which 
must avoid e\ will prefer to return to the origin from — ei, so the constraint is greater when 
this possibility is disallowed. In fact this intuition is borne out by the numerical results, and 
as can be seen from Table 10 it is also evident from the 1/d expansion. From Table 10, we 
have 

9 rl — 1 

(2jl) = m = _ _ _ ^ _ _ i24 ^ 4 

C 2 {ei} -(0,0;^_) 

which has the correct coefficient of order s 2 , and as a bonus only misses the correct coefficient 
of order s 3 by 2. It is also clear from Table 10 that the expansion for fji 2 ' 2 ^ is identical to 
that of ^ 2,1 \ through the order shown in (6.24). 

The optimized bounds of Section 2.5 do not improve substantially on (6.24): we find that 
for the optimized bound using k = 2 and r = 2 the terms up to and including order s 3 are 
as in (6.24), while the coefficient of s 4 is improved slightly from —124 to —122. 

In general we do not expect that there will be improvements below order s 5 when k is 
increased beyond 1, for a given memory. For example, for k = 2 the number of hexagons 
which pass through a specific next-nearest neighbour of the origin is only O(d) , and hence 
when multiplied by f3 6 is an order s 5 effect. This can be viewed as a partial explanation of 
the small size of the improvements observed in the numerical bounds as \A\ increases beyond 
1. 



6.4 1/d expansion for the Ising-model lower bounds 

Let us first look at the bound (5.5) obtained by using comparison to the Ising model together 
with the infrared bound on J c ,ising- This improves our basic r = 0, k = bound by virtue 
of the coth, and yields 

// > s" 1 - 1 - §s - f s 2 - ±§V + 0(s 4 ) . (6.25) 

We see that the coth makes only a slight improvement in the term of order s, and does not 
achieve the correct coefficient —1. 

The method based on comparison to the Ising model is potentially more powerful than 
this, if one could get a better bound on J c ,ising- The best possible result is obtained by 
inserting the exact (nonrigorous) expansion 

Jcjsing = t -i 4 13 ?, 979 3 2009 4 1 (6.26) 

* 1 3 d 3 * 45 ^ 15 * T 

derived by Fisher and Gaunt [26] into (5.3); the result is 

// > coth(J Cj/sm3 ) = s' 1 - 1 - s -4s 2 - 21s 3 - ^s 4 + ■ ■ ■ . (6.27) 
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This achieves the correct first three terms, but does not do as well as (6.24) on the term 
of order s 2 . Of course, it is a highly nontrivial open problem to prove rigorously such good 
bounds on J c j stng . 



A Simple random walk 

A.l Numerical evaluation of the two-point function 

To evaluate our lower bounds on // numerically, it is necessary to know the numerical values 
of the critical simple-random-walk two-point function Co(0, x; ^) for various values of x. An 
effective method of evaluating this quantity (as well as the subcritical two-point function) 
numerically to high precision, with rigorous error bounds, has been described in considerable 
detail in Appendix B of [4]. We now summarize briefly how the calculation goes in the critical 
case. 

We begin with the integral formula 

/d d k ^/ik-x 
7 w — • (A.l) 



where 



We then apply the identity 



_ 1 d 

D(k) = -J2 cos h. (A. 2) 

8 = 1 



1 ,' " e~ At dt , (A. 3) 



which is valid for A > 0, with A = 1 — D(k). This factorizes the integrals over £4, . . . , kd, 
which can then be performed to give 

00 d 

C o (0,x;^) = dj f[f\ x S)dt, (A.4) 

8=1 

where 

f N (z) = e- z I N (z) = — f e - z(1 - cose ^ cos NO dO (A.5) 

— 7T 

and In(z) is the modified Bessel function. This transforms the o?- dimensional integral (A.l) 
into a one-dimensional integral over a semi-infinite interval. 

The integrand in (A.4) decays as t — > 00 only as a power. We speed up the decay by 
making the change of variables t = e", obtaining 

00 

C o (0,x;^) = d [ F(u)du (A. 6) 
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A n (x) 


Co(0,x;±) 


X 


d = 2 


a = 3 


a = 4 


a = 5 


a = 6 


( C\\ 


n 

u 


1.5163860744 


1.2394671218 


1.1563081248 


1.1169633732 


(\ \ 
\ l ) 


1 

1 


0.5163860744 


0.2394671218 


0.1563081248 


0.1169633732 




1 973930^4-/17 

1 .Z 1 OZO^O^l^l 1 


0.3311486174 


0.1017176302 


0.0474085960 


0.0271774706 


\ L i L i L ) 




0.2614701416 


0.0618723811 


0.0222517907 


0.0100651251 


i^ ^ ^ X\ 
I 1 , 1 , 1 , 1 ,! 







0.0447274307 


0.0133523237 


0.0049990827 


^1,1,1,1,1; 







— 


0.0092253734 


0.0029872751 


\ z ) 


i .^ooozuy iuo 


0.2573359025 


0.0659640719 


0.0275043553 


0.0148223998 


(9 1 


1 ^/Ifizl^QnSQ^ 

1 .rJ^lU^l i VVJOVO 


0.2155896361 


0.0436586366 


0.0139794831 


0.0058409498 






0.1917916659 


0.0334570990 


0.0089609415 


0.0030848645 


(9 1 1 1 

[Z, 1,1, 1J 







0.0275824802 


0.0065163319 


0.0019448478 


(9 9\ 
\ Zj i Zj ) 


1 fiQ7fi^979fi3 
l.Uj 1 UOZ 1 zuo 


0.1683310508 


0.0259898362 


0.0062387819 


0.0019599801 


(9 9 1 




0.1569524280 


0.0221867673 


0.0047601264 


0.0013047592 




1 791 1 9^4639 


0.1652707962 


0.0262936339 


0.0068995628 


0.0024959268 


(3,1) 


1.7615031763 


0.1531389140 


0.0217691587 


0.0048774493 


0.0014526312 


(3,1,1) 




0.1441957255 


0.0189286425 


0.0038130772 


0.0009927439 


(3,2) 


1.8488263632 


0.1324510884 


0.0159271735 


0.0029340473 


0.0006998942 


(4) 


1.9079745896 


0.1217332189 


0.0137700477 


0.0024716782 


0.0006024105 


(4,1) 


1.9295817894 


0.1171305125 


0.0125592552 


0.0020829365 


0.0004528519 


(5) 


2.0516093163 


0.0966064672 


0.0085112166 


0.0011537265 


0.0002044795 



Table 4: Numerical values of A (x) for d = 2, and of Co(0,x; ^) for d = 3,4,5,6, for the 
values of x needed to compute the lower bounds on The values are rounded to ten digits 
after the decimal point. Final components of x which are not shown are equal to zero. 



where 

F(u) = e u Hf lxtl (e u ). (A.7) 

8 = 1 

We then use standard methods to bound the difference between the integral (A. 6) and an 
infinite Riemann sum, truncate the Riemann sum with bounds on the omitted tails, and 
evaluate the resulting finite Riemann sum, using the large- z asymptotic expansion for the 
modified Bessel function to deal with large t and a truncated Taylor series for the modified 
Bessel function to deal with the remaining t. We also take into account possible round-off 
errors in the numerical computations. The result is that we obtain the values in Table 4. 
The values of A (x) for d = 2, which are known exactly, have been computed using the 
algorithm described in Section 15 of [15]. 

Table 5 gives numerical values of C^(0,0;^), computed as described in Section 3.1. 
Table 6 gives the values of C<f (0, 0; ^777), computed using the method of Section 3.2. Table 7 
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A 


J = 2 


d = 3 


d = 4 


d = 5 


d = 6 


(PS 


OO 


1 .OIUOOU 


1 93Q4R7 


1 . 10UOUO 


1 . 1 iuyuo 


{p \ 
i e l/ 


o 
z 


1 .0'lUOO 1 


1 1 Q3909 
i . lyozuz 


1 1 3^1 7Q 


1 1 0471 ^ 

1 . lU^ 1 10 


if- 9p- 1 
t e l, ze lj 


1 .OrJUUOU 


1 33^84 
1 .ooooo^t 




1 13^143 


1 1 0470Q 




1 7^Q10 




1 19069^ 


1 1^4^70 


1 104^14 


{ei,2ei,3ei} 


1.789955 


1.333592 


1.192796 


1.135138 


1.104709 


{ei,2ei,2ei + e 2 } 


1.767994 


1.330950 


1.192528 


1.135103 


1.104703 


{ei,ei + e 2 ,2ei + e 2 } 


1.695733 


1.320170 


1.190465 


1.134555 


1.104512 


{ei, ei + e 2 , e x + 2e 2 } 


1.658351 


1.317957 


1.190277 


1.134530 


1.104507 


{ei, ei + e 2 , e x + e 2 + e 3 } 




1.314688 


1.189764 


1.134435 


1.104485 


{ei, e x + e 2 , e 2 } 


1.357421 


1.226256 


1.153726 


1.115672 


1.093040 



Table 5: Numerical values of Qf(0,0; ^) for d = 2,3,4,5,6, for various choices of A. The 
values are rounded to 6 digits after the decimal point. 



gives numerical values for (7^' e (0, 0; ^777), computed using the method of Section 3.3. 

A. 2 \jd expansions 

This section contains 1/d expansions for several relevant quantities. As before, we use the 
shorthand s = First, the following are obtained by directly integrating powers of cosines. 



d d k 
(2ttY 



D(ky 



s 

3s 2 - 3s 3 

15s 3 - 45s 4 + 40s 5 
105s 4 - 630s 5 + 1435s 6 
945s 5 + 0(s 6 ) 
0(s 6 



1155s' 



2) 
4) 
6) 
8) 
10) 



d d k 
(27r) a 



D(k) n cos 2^ 



d d k 
(2ttY 



D(k) n cos k\ cos k 2 





(n > 12 




(n = 2) 


6s 3 - 8s 4 


(n = 4) 


< 45s 4 - 165s 5 + 165s 6 


(n = 6) 


420s 5 + 0(s 6 ) 


(n = 8) 


. 0(s 6 ) 


(n > 10) 


' 2s 2 


(n = 2) 


12s 3 - 24s 4 


(n = 4) 


< 90s 4 - 450s 5 + 660s 6 


(n = 6) 


840s 5 + 0(s 6 ) 


(n = 8) 


. 0(s 6 ) 


(n > 10) 



(A.8) 



(A.9) 



(A.10) 
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A 


J = 2 


d = 3 


d = 4 


d = 5 


d = 6 





(DO 


1.213109 


1.062400 


1.027829 


1.015421 


M 


1.751695 


1.125805 


1.045320 


1.021995 


1.012778 


{ei,2ei} 


1.564686 


1.121796 


1.045056 


1.021965 


1.012773 


{ei, e x + e 2 } 


1.481256 


1.113713 


1.043739 


1.021661 


1.012681 


{ei,2ei,3ei} 


1.485308 


1.120113 


1.044989 


1.021961 


1.012773 


{ei,2ei,2ei + e 2 } 


1.465260 


1.118079 


1.044790 


1.021935 


1.012769 


{ei,ei + e 2 ,2ei + e 2 } 


1.423896 


1.111657 


1.043606 


1.021649 


1.012680 


{ei, ei + e 2 , e x + 2e 2 } 


1.388210 


1.110045 


1.043471 


1.021631 


1.012677 


{ei, e x + e 2 , e x + e 2 + e 3 } 




1.107650 


1.043096 


1.021561 


1.012660 


{ei, ei + e 2 , e 2 } 


1.211771 


1.071615 


1.031850 


1.017013 


1.010433 



Table 6: Numerical values of C^(0, 0; -^jzi) f° r = 2, 3, 4, 5, 6, for various choices of A. The 
values are rounded to 6 digits after the decimal point. 



A, e 


J = 2 


d = 3 


d = 4 


d = 5 


d = 6 


0,e 


oo 


1.177591 


1.054600 


1.025047 


1.014136 


{ei}, -ei 


1.304467 


1.097370 


1.038115 


1.019313 


1.011519 


{ei},e 2 


1.349274 


1.101462 


1.038967 


1.019581 


1.011627 


{ei, 2ei}, -ei 


1.193791 


1.093963 


1.037871 


1.019285 


1.011514 


{ei, 2ei}, e 2 


1.238380 


1.098305 


1.038744 


1.019554 


1.011622 


{ei, ei + e 2 }, -e t 


1.140666 


1.086889 


1.036639 


1.018992 


1.011424 


{ei, ei + e 2 }, e 2 


1.210882 


1.094439 


1.038093 


1.019408 


1.011579 


{ei, ei + e 2 }, -e 2 


1.170686 


1.090659 


1.037469 


1.019256 


1.011531 


{ei, ei + e 2 }, e 3 




1.091432 


1.037558 


1.019272 


1.011535 



Table 7: Numerical values of C 2 ' e (0, 0; -^—[) f° r ^ = 2,3, 4, 5, 6, for various choices of A, e. 
The values are rounded to 6 digits after the decimal point. 
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31 



A 

/ 1 


U l U 5 U > 2 J 




1 + 5 + 3s + 12s + 60s + 355s + C(s J 


{ei} 


1 i io2iT3ior'4ioir'5i f~\ ( 6 \ 

1 + 5 + 2s + 7s + 35s + 215s + U(s ) 


{e 1 ,2e 1 \ 


1 i io2iT3ior'4ioir'5i f~\ ( 6 \ 

1 + s + 2s + 7s + oos + 215s + U(s ) 


r i i 

{ei, ei + e 2 ) 


1 i io2iT3iO/i4ior\o5i f~\ ( 6 \ 

1 + s + 2s + 7s + 34s + 202s + U(s ) 


fool 


1 i io2iT3ior'4ioir'5i f~\ ( 6 \ 

1 + s + 2s + 7s + 35s + 215s + U(s ) 


{ei, 2ei, 2ei + e 2 j 


1 + s + zs + is + o5s + zl5s + U(s ) 


{ei,ei + e 2 ,2ei + e 2 } 


1 + s + 2s 2 + 7s 3 + 34s 4 + 202s 5 + 0(s 6 ) 


{ei, ei + e 2 , e x + 2e 2 } 


1 + s + 2s 2 + 7s 3 + 34s 4 + 202s 5 + 0(s 6 ) 


{ei, ei + e 2 , e x + e 2 + e 3 } 


1 + s + 2s 2 + 7s 3 + 34s 4 + 202s 5 + 0(s 6 ) 


{ei, e x + e 2 , e 2 } 


1 + s + s 2 + 2s 3 + 14s 4 + 115s 5 + 0(s 6 ) 



Table 8: 1/d expansions for C^(0,0; ^) for various choices of A. 



Finally, we turn to the 1/d expansions for loop-generating functions with taboo set. 
Beginning with the 1/d expansions for Co(y, x; ^) = Iw(x — y) given above, and then using 
the recursion (3.1), we obtain the 1/d expansions for C^(0,0; ^) given in Table 8. The 1/d 
expansions for memory-2 quantities are more difficult to handle. As described under (3.13), 
we first solve (3.13) for {C^^iy , b + /)}|/|=i, for fixed y. At this stage, special care is 
needed to make full use of symmetry, since naively (3.13) is a system of equations for 2d 
unknowns (and here d /* oo). By symmetry, we can reduce (3.13) to a system of equations 
for a number of unknowns which is uniformly bounded in d (at least for small \A\ and |j/|, \x\), 
and obtain the results in Table 9. Then, using (3.18), we compute the 1/d expansions for 
C^' e (0,0) given in Table 10. 
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A 


C 2 A (0,0;^t) 





1 + s 2 + 7s 3 + 43s 4 + 278s 5 + 0(s 6 ) 


M 


1 + s 2 + 5s 3 + 29s 4 + 188s 5 + 0(s 6 ) 


{ei,2ei} 


1 + s 2 + 5s 3 + 29s 4 + 188s 5 + 0(s 6 ) 


{ei, e x + e 2 } 


1 + s 2 + 5s 3 + 29s 4 + 182s 5 + 0(s 6 ) 



Table 9: \jd expansions for C^(0,0; j^r[) for various choices of A. 



A, e 


^(0,0;^) 


0,e 


1 + s 2 + 6s 3 + 36s 4 + 235s 5 + 0(s 6 ) 


{ei}, -ei 


1 + s 2 + 4s 3 + 22s 4 + 146s 5 + 0(s 6 ) 


{ei},e 2 


1 + s 2 + 4s 3 + 23s 4 + 154s 5 + 0(s 6 ) 


{ei, 2ei}, -ei 


1 + s 2 + 4s 3 + 22s 4 + 146s 5 + 0(s 6 ) 


{ei, 2ei}, e 2 


1 + s 2 + 4s 3 + 23s 4 + 154s 5 + 0(s 6 ) 


{ei, ei + e 2 }, -e\ 


1 + s 2 + 4s 3 + 22s 4 + 140s 5 + 0(s 6 ) 


{ei, ei + e 2 }, e 2 


1 + s 2 + 4s 3 + 23s 4 + 152s 5 + 0(s 6 ) 


{ei, ei + e 2 }, -e 2 


1 + s 2 + 4s 3 + 23s 4 + 148s 5 + 0(s 6 ) 


{ei, ei + e 2 }, e 3 


1 + s 2 + 4s 3 + 23s 4 + 148s 5 + 0(s 6 ) 



Table 10: 1/d expansions for C 2 ' e (0, 0; -^—[) f° r various choices of A, e. 
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